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Abstract 

Spin-wave instability is studied analytically in the case of parallel pumping for thin 
films under external field perpendicular to the film plane. It is examined whether 
the instability threshold derived from only the Landau-Lifshitz (LL) equation can 
explain experimental instability threshold without using the microscopically-derived 
spin- wave line width AH^, which is conventionally used. It is revealed that the but- 
terfly curve cannot be explained from only the LL equation at least in an analytical 
way. By contrast, for the case of perpendicular pumping, the Suhl instability was 
well explained from the LL equation. The difference between the two cases comes 
from the nonlinear terms describing the relaxation of spin waves. It is suggested 
how the nonlinear terms in the LL equation should be related to AHk for parallel 
pumping. 
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1 Introduction 



Spin-wave instability is a nonlinear phenomena which is dominated by a bal- 
ance between excitation and relaxation of spin waves. In early 1950s, anoma- 
lous phenomena incompatible with the conventional theory were observed 
in perpendicular-pumping ferromagnetic resonance (FMR) experiments [1,2]. 
The phenomena was successfully explained by Suhl's theory [3]. In the theory, 
the coupling between the uniform mode {k = 0) and other spin-wave modes 
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(k 7^ 0) plays a crucial role. In perpendicular pumping, the microwave mag- 
netic field excites the uniform mode. Energy is fed from a A; = mode or a 
pair of A; = modes to a pair of spin waves k and —k {k ^ 0). The pair of 
spin waves which is most strongly coupled to the uniform mode grows expo- 
nentially when the intensity of the pumping field exceeds a certain threshold. 
Such processes can be well described by using only the Landau-Lifshitz (LL) 
equation governing the magnetization dynamics. 

On the other hand, the spin-wave instability for parallel pumping has been 
investigated in order to study relaxation phenomena of spin waves with def- 
inite wave vectors since the first experimental studies in 1960s [4,5]. In the 
parallel pumping experiments, the microwave magnetic field is applied par- 
allel to the equilibrium magnetization. Therefore the longitudinal component 
of a standing spin-wave, which consists of a pair of traveling spin-waves k 
and —k, interacts with the microwave magnetic field. When the power of the 
pumping field exceeds a certain threshold, the spin waves k and —k grow 
exponentially. The threshold depends on the balance between excitation and 
relaxation of spin waves. The relaxation of spin waves can be described by the 
spin-wave line width AHi^ which is derived from microscopic magnon-magnon 
scattering theories. To explain the experimentally-observed spin-wave instabil- 
ity threshold, i.e. so-called butterfiy curve, AHk is used as a factor of damping 
in most of conventional theories. For example, Patton et al. proposed some 
trial AHk functions and successfully explained butterfiy curves of the insta- 
bility threshold [9]. In contrast to the Suhl instability in the perpendicular 
pumping, however, there has been no attempt to use only the LL equation for 
explaining the experimental butterfly curve. 

In the 1980s, there occurred a renaissance on studies of FMR [6-8]. High- 
resolution experiments were carried out for spin-wave nonlinear dynamics in 
yttrium iron garnet (YIG) fllms and spheres in parallel and perpendicular 
pumping methods. For certain high power beyond the Suhl instability thresh- 
old, the interaction between excited spin-wave modes lead to various dynam- 
ical phenomena such as chaos and turbulence. Since then the studies on spin 
wave instabilities have acquired a renewed interest up to now [10-12]. Despite 
these pioneering experiments, the corresponding theoretical investigations for 
parallel pumping leave several crucial questions. In particular, no dynamical 
equation leading to auto-oscillations and chaos [13] and pattern formations [14] 
can touch on even a qualitative feature of the experimental butterfiy curve for 
the instability threshold. In the case of parallel pumping, therefore, it is im- 
portant to analyze the butterfly curve from a viewpoint of modern nonlinear 
dynamics. 

In this paper, we try to derive analytically the spin-wave instability threshold 
in ferromagnetic thin films for parallel pumping without using the spin-wave 
line width AiJ^ which comes from microscopic theories. By using only the 
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LL equation, besides the linear instability analysis, the nonlinear instability 
is examined for a simple case. It will be shown that the instability threshold 
obtained in that way cannot explain experimental butterfly curves. We will 
reveal why our nonlinear instability analysis fails and suggest what terms 
arising from the LL equation should correspond to Aif^. 



2 Equations of motion for the spin-wave amplitudes 



The LL equation is given by 

-dtM ^MxH^ii- -^M X (M X iJeff). (1) 

7 J^o 



Here, M{r) is a magnetization field and Mq is its value in thermal equilib- 
rium; 7 is the gyromagnetic ratio (7 < for spins); A is a damping constant 
parameter (A < in this case); Jfeff is the effective magnetic field: 

H^s = DV^M + H'^ + H'' + Ho + hcosujt. (2) 



The first term comes from the exchange interaction; the second term is the 
demagnetization term; the third term is the anisotropy field, which is omitted 
below for convenience; the fourth and fifth terms are the external static and 
pumping fields, respectively. We assume that Hq/ /h and these fields are par- 
allel to the z axis. The demagnetization field is given as the gradient of the 
magneto-static potential (p: 

= -V0. (3) 



The magneto-static potential obeys the Poisson equation: 
47rV • M, inside the sample; 



(4) 

0, outside the sample. 



Throughout the text, we call a set of Eqs. (l)-(4) as the LL equation. Now we 
introduce the normalized magnetization S = M/Mq, considering the length 
of M is invariant: Mq = \M\. It is convenient to project S stereographically 
onto a complex variable ip{r, t) [15]: 

Sx + iSy 
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where 



(6) 



In terms of ip, Eq. (1) is rewritten as 



dtip - i(l - iA) 7 <^ DMo 



1 + ipijj* 

- ^{1 + ip*^)dy(f) + {d^cf)- Ho- h cos (jjt)^^ = 0. 
Inside the sample, (p satisfies 



(7) 



AttMq 



1.2 



+ i[(l + r)dyij* - (1 + ij*')dyij] - 2{i)d,^* + ^*d, 



(8) 



The demagnetizing field is affected by boundary conditions. We consider a 
film of thickness d infinitely extended in the x-y plane under the external filed 
applied parallel to z axis. We assume unpinned surface spins, which satisfy 
Neumann-like boundary conditions. 



(9) 



namely. 



\z=±d/2 



0. 



(10) 



Introducing the dimensionless time and space units [14], 
t 



t 



47r|7|Mc 



■, r — > rd, 



we obtain the linearized equations of motion corresponding to Eqs. (7) and 
(8): 



dtiP + i{l-iX) 



I'V'ijj -^{dx + idy)^ + {d,^ -UH-uJh cos ujpt)ilj 



0,(12) 



= {d^ - idy)ilj + {d^ + idy)ilj* when -| < z < |, 



(13) 
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where 



4^^' "^"'Mfo' "^""'Mfo' '^^"i^^MoH 



First of all, we consider the undriven case (i.e. cuh = 0) and expand ip{r,t) 
and ip*{r, t) so that they fulfills the boundary conditions (10). For even modes, 
kz — 2m7r (m: integer). 



ijj{r, t)=J2 afc(t) e'(*^-^+*^-2') cos k^z 

k 

jj;*(r, t) = Y^ a-k{t) e'^'^-^+'^yy^ cos k^z. (15) 
fc 

For odd modes, k^ — (2m + l)7r (m: integer), 

k 

j/j* (r, i) = - ^ alfe(i) e'('=-^+'=^^) sin (16) 
fe 

Using the expansions (15) and (16), we obtain solutions of Eq. (13): 



Qi{kxX+kyy) 

$(r) = — i X] 72 (^-^fe + k+(^*-k)fk cos kzZ for even modes; 

fc ^ 

Qi{kxX+kyy) 

$(r)==— iV" — {k-Gk — k+a'iujfkSmkzZ for odd modes. (17) 

V ^ 

Here k+ — kx + iky, k- — kx — iky, k± — ■\Jkl + ky and 



/fe=l-(l-e-'=-)|^ when A;, -0; 
/fe = l-2(l-e-'=-)|^ whenA;,^0. 

(18) 

The detailed derivation of these solutions is shown in Appendix A and Ref . [16] . 

The derivatives 9^$ and dy^ are calculated from Eq. (17). However, 9^$ = 
—H^/AttMo is approximated with the value of uniform magnetization, k = 0: 
= —AttN^Mz, where is a demagnetizing factor. In this case, = 1. 
Considering ~ Mq, we have d^^ — 1. Then Eq. (12) for = is rewritten 
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as 



dti^ + i(l - iA) 



rV^^ _ + idy)^ + (1 - ujh)iI^ 



0. 



(19) 



From Eqs. (15)-(19), the linearized equation of motion of is found to be 
dtak + i(l - iA)Afeafc + i(l - \\)Buaij^ = 0, (20) 

where 



1 

2* 



Ak = 


1 - l^k 




(-1)' 


Bk — 




2 


sin 9 k = 







e^'^'^/feSin^^fe, 



Here, n is an integer and kz = Tm. 



tl. _|_ { 

k\ k\ 



[21] 



Equation (20) represents two coupled harmonic oscillators afc, a*_^. and can be 
diagonalized by means of the Holstein-Primakoff transformation: 



flfc — T^kbk — fJ'kb*_k 

a*_k = i^kb-k - I^U>k, 



(22) 



where 



i/fe = cosh ^ , /ifc = e^^'= sinh ^ , 



Xfe 



cosh Xfc ■■ 



I -4 J 



1/2- 



A|-(1 + A2)|i?fc|2_ 
After the transformation (22), Eq. (20) becomes 
dfbk - i{uJk + ir]k)bk = 0, 



(23) 



(24) 



where 



cvl^Al-{l + X')\Bk\\ 
rik = A^fe 



(25) 
(26) 



Equations (25) and (26) express the dispersion relation and a damping rate 
for the spin wave, respectively. 
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3 Linear instability 



We will first analyze the linear instability for the spin wave under the driving 
field. When the pumping field h cos cut is applied, the equation of motion of 
Ofe corresponding to Eq. (20) becomes 

dfttk + i(l - iX)Akak + i(l - iX)Bka*_^ - i(l - iX)uJhCos{upt)ak = 0. (27) 
After the the transformation (22), we substitute the following equations into 



since the resonance occurs at ujk = When only terms contributing to 

the resonance are left, the variable 6^ satisfies 



Eq. (27), 



bk{t)^bl{t)exp[i{u;j2)t-rikt], 
b*_^{t) = b%{t) exp[-i(a;p/2)i - rikt], 



(28) 




(29) 



where 



Pk\ = UhV 1 + A2 



(30) 



2uk 



Therefore, the exponentially increasing solution for bk. is 




(31) 



where 



Pk\ > Vk- 



(32) 



From Eqs. (30) and (31), the instability threshold oo'fl'^'^ is given as 




(33) 



Let us rewrite Eq. (33) with use of Eqs. (25) and (26): 



By using Eqs. (25) and (34), we might sec a theoretical butterfly curve, which 
is an instabihty threshold curve plotted against the static field. However, the 
butterfly curve calculated in that way is found to be totally different from 
experimental ones. Typical butterfly curves observed in experiments have a 
cusp at a certain static held: as the static held increasing, a;"'* flrst decreases 
up to the cusp point and then increases above that point. Here let us focus 
on the case for static fields below the cusp point. For those static field, the 
minimum threshold mode corresponds to a spin wave propagating with 9^ — 
7r/2, i.e. parallel to the film plane. As the static field increases, the wave vector 
k of the threshold modes decreases, and A; ~ at the cusp. On the contrary, 
from Eq. (21), one notices that cj™* in Eq. (34) depends only on fk when 9k 
is fixed and that it grows as the static field increases up to the cusp point. In 
fact, for 9k = 7r/2, 

f^^l-(l-e-')^ (35) 

because kz = and k± = k. As k ^ 0, fk and then ujf^^^ oo. Such 
a divergence of uu^^^^ — >■ oo has not been observed in experiments. Namely, 
the theoretical threshold u;^^^ in Eq. (34) fails to explain the experimental 
butterfly curve. 

One may say that the linear analysis is not sufficient to discuss the spin-wave 
instability threshold. In the next section, several nonlinear terms are included 
in equation of motion, and a nonlinear instability is analyzed with expectation 
to overcome the above problem. 



4 Nonlinear instability 

Let us rewrite Eqs. (7) and (8) up to third order of ip. Then we have 

dtijj + i(l - iA) [V'^ - 2rWf] - + i9,)$ + l^Jj^d:, - idy)^ 

+ ((9^$ -uiH -uih COS u)pt)tp} , (36) 

and 

Our interest lies in the region of static flelds below the cusp point, where 
9k — 7^/2. Then the expansion of t/j is 
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k 

ijj* (r, t) = J2 ci-kit) e'(^-"+'=«s') . (38) 
fe 

Substituting Eq. (38) into Eq. (37), we have 
fe 

felfe2fe3 

X [(A;3+ + 2ki+)a*_f^^a*_f^^ak3 + (ks- + 2k2-)akiak^a*_^.^ . (39) 

It is convenient to restrict wave vectors for solving this equation. When ki + 
^2 + ^3 = k, the solution of Eq. (39) is 



Q\(kxX+kyy) 

k 

+ E [(^3+ + 2kl+)a*_^.y_^^^ak^ + (/c3_ + 2k2-)ak^aky_k^ } /fe, (40) 

felfe2fe3 

where 

1 _ p-fcx 

/fe = 1 - (41) 



The derivation of these equations is similar to that of the linearized equations. 

The nonlinear terms contributing to the resonance have spin waves whose wave 
vector is fc or — fc. When the linear terms have fc, a possible combination of 
fci, fc2 and fca is (fc, fc', — fc'). Here, however, we assume fci, fc2, fcs = ±fc. This 
assumption, which forbids the multi-mode couplings of spin waves, incorpo- 
rates the substantial nonhnear terms. Then, from Eqs. (36), (38) and (40), we 
obtain 



dtttk + i(l - iA) jAfeafe -Bfeal^, - Uh cos{uj^t)ak 

+ 2(Cfc fk)aka_kal.k + ("Cfe + /fe)«fe«fe«fc 

+ i/fee-'''^fe akaka^k + i/fee''<^fc(2alfca;^afc + alf,a*_^a.k)} = W 

where 



Ak^l- fk'^ -ujh- ^/fe, 
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(43) 



As in the case of linear instability, we perform the Holstein-Primakoff trans- 
formation (22) and substitute Eq. (28). When only terms contributing to the 
resonance are left, the equations of motion for and 6°*^, are 

dtbl -i(u;k-^+ ink) hi + \i,:\hl\X + 2Kk\hW%l - pkb'>_\ = 0, 

9,6- +i K - ^ - ivk) b'L*^ - iaKkl'b'^-k - '^KlKl'b'^-k - PIK = m) 

where 



ik—Ck 



_ 1 _ ^ 



\ 3|i?fep(l + A^) ir/fe 
[ 2a;fc a;fc 



3|Afc|/fcA2 3iA/fc' 



4u;i 



Ck — Ck 



+ 1 + ^ 



/ 3|EfcP(l + A^) i?7fc 3| Afcl/fcA^ 3iA/fe ' 



'k 

The detailed derivation is given in Appendix B. 



4u;t 



(45) 



The time- independent solutions (fixed points) of Eq. (44) ^b% — W and b°*^ — 
W , are given by 



W 



2 _ r;felm(a + 2^) ± [vUMCk + 2a)}' + 1^ + 2^1' {\Pkf - vl) 



1/2 



la + 2Cfe| 



-.(46) 



Since the inside of the square root in Eq. (46) should be positive. 



, ,2 . [Re(a+2a)F , 

l^fel — 17 , ^> |2 'Ik- 



ia+2Cfei 



(47) 



W 



> 0. 



When Ipfcp > r]1, the double sign changes into the plus sign because 

By using b1 = W + bi and b"*^ = W* + 62 in Eq. (44) , the hnear equations of 
bi and 62 are written as 




Vk + 2i(a + 2Cfe) 



W 



KCk + 2Ck)w-pk 



Ka + ^Ck)w*' -Pi Vk- 2i(a + 2a) 



w 



hi 



,(48) 
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where we assumed 6* = 62 and 63 — The eigenvalues of the matrix are 
given as 



A+ = -rjk + 2Im(efc + 2a) 
A_ = -rjk + 2Im(a + 2Cfe) 



W 



+ 



w 



r/2-4{Re(a + 2a)F 
ry2-4{Re(a + 2a)r 



W 



1/2 



W 



nl/2 



■ (49) 



Since A+ > A_, the instabihty condition is A+ > 0. However, this condition 
does not lead to any nontrivial condition. Then the actual instability condition 
is Eq. (47). The instability threshold now becomes 



crit 



ujpr]k |Re(a + 2Cfc)| 

\Bk\ VTT>? la + 2Cfe| 



(50) 



However, this equation cannot explain the experimental butterfly curve of the 
instability threshold, either. The factor |Re(a + 2Cfc)|/|a + 2Cfc| is almost 
constant for the region where our calculation is performed, and we cannot 
overcome the problem of divergence cuf^^^ — > 00 as A; — > 0. 



5 Discussion 



We have tried to derive the spin- wave instability threshold for parallel pumping 
from the LL equation. Neither the linear nor nonlinear analysis has proved 
to explain the feature of experimental butterfly curves of the threshold. By 
contrast, the Suhl instability for perpendicular pumping can be well explained 
from the LL equation [3] . The successful explanation of Suhl instabihty is based 
on the fact that the couplings between the uniform mode and other spin-wave 
modes are dominant for perpendicular pumping. In parallel pumping, however, 
the uniform mode is not excited. That is the reason why the Suhl instability 
cannot be applied to the present case. 

In order to explain the instability threshold for parallel pumping, it is necessary 
to deal with the spin-wave relaxation in a proper way. Conventionally, the spin- 
wave line width AH^, which is calculated from the microscopic theory, is used 
to describe the spin- wave relaxation. In Ref. [9], the /c-dependence of Aifj. 
plays a vital role in obtaining the butterfly curve. Such a /c-dependence of 
AHk should correspond to the multi-mode spin-wave coupling, which is not 
considered in the present treatment of the LL equation. However, the equation 
of motion with multi-mode couplings is too complex to be solved analytically. 
That may be the reason why the instability threshold for parallel pumping 
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has never been derived in the context of nonhnear dynamics by using only the 
LL equation. 
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A Demagnetization field 

Substituting Eqs. (15) and (16) into the Poisson equation, we have 



V2$ = i^e'('=-^+'=^2')(A;_afc - k+a*_^^) sink.z for odd modes. (A.l) 



= i e'('''^^+''«^)(A;_afc + k+a*_J cos k^z for even modes; 



k 



k 



Generally, for the Poisson equation given as 



V^(/)(r) = -47rp(r), 



(A.2) 



the solution is 




(A.3) 



By using Eq. (A.3), the solutions of Eq. (A.l) are written as 




for even mode; 



k 




sin kzZ' 



for odd mode. 



(A.4) 



After integration, Eq. (A.4) becomes 
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$(r) = -iYl U (^-«fc + ^+«-fe) cos k^z - (-1)"" e"''-^/^ cosh fc^^ 

for even mode; 

pi{ka:X+kyy) 

$(r) = -iV (k^uk - k+a*_u) sin k,z - (-1)'^ e'''^/'^ sinhk^z 

k ^ L J 

for odd mode. (A. 5) 

Now let us introduce the projection F{z) of an arbitrary function f{z) onto 
cos kzz: 



1/2 1/2 

F{z) — j dzf{z) COS kzZ / J dzcos'^kzZ, (A. 6) 

-1/2 ' -1/2 



for —1/2 < z < 1/2. The projection onto sin k^z is also defined in the same 
way. Applying the above projection, we obtain 

Qi(kxX+kyy) 

$(r) = — iT^ — {k-ttk + k+a*_u)fk cos kzZ for even modes; 



k 



^i{kxX+kyy) 

$(r) = — iT^ — (k_ak — k+a*_u)fkSinkzZ for odd modes. (A. 7) 

k^ 



Here, 



/fe = l-(l-e-'=^)^ whenA;, = 0; 

/fe = l-2(l-e-'=-)|^ when A;, ^0. (A.8) 



B Nonlinecir equation of motion for b% 



First of all, let the terms of Eq. (42) be represented by Pi,fc, P2,fe, -Ps.fe, P4,k 
and P5,fe: 



Pi,k = dtttk + i(l - iA) (^AfcOfc + ^fcCilfc - i^h cos{upt)akj 
P2,k = 2i(l - iA) (Cfc + fk)aka^ka*_k 
-P3,fe = i(l - iA)(-Cfc + fk)akakal 
PA,k = |i(l - iA)/fce-2i<^fe akttka-k 

P,,k = ^i{l - iA)/fce2'^fc(2a*_fca;^afc + a*_fca*_fca_fc). (B.l) 
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Next, we introduce a new quantity Ri^k = i^kPi,k + l^kPi-k for z = 1, . . . , 5. 
After the Holstein-PrimakofT transformation (22), we use Eq. (28) and remove 
the terms which do not contribute to the resonance. Then, we have 



where 



pk = i(l - iA)a;^- — . (B.3) 
Zojk 



The terms in Eq. (B.2) satisfy the following: 

Rl,k + R2,k + -Rs.fc + -^4,*; + -Rs.fc = 0. (B.4) 

Equation (44) is obtained from Eqs. (B.2) and (B.4). 
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